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We obtained numerical and closed-form analytic expressions for finite-temperature plasmon dis¬ 
persion relations for intrinsic graphene in the presence of a finite energy gap in the energy spectrum. 
The calculations were carried out using the random-phase approximation. The analytic results have 
been derived in the high temperature regime and long-wavelength limit. We have found that the 
plasmon damping rate decreases in the presence of a band gap. Our method of calculation could 
also be applied to silicene and other buckled honeycomb lattice structures. The finite-temperature 
plasmon dispersion relations are presented when a single graphene layer is Coulomb coupled to a 
semi inhnite conductor. Both cases of gapless and gapped monolayer graphene have been inves¬ 
tigated when a thick substrate is in their proximity. Both the plasmon excitation frequency and 
damping rate are linear functions of the in-plane wave vector in the long wavelength limit when a 
monolayer interacts with a conducting substrate which is not the case for free-standing pristine or 
gapped graphene. 


I. INTRODUCTION 

Plasmons in grp ahene represent one of the most interesting and actively studied fields, from both fundamentally 
theoreticaP^^^EHin] and experimental or technological points of vievE®. Graphene plasmons are especially important, 
partially because of their versatile frequency range which may be adjusted by varying the doping concentrations of 
its mobile carriers and ener gy b and gap. Gonsequently, graphene has several potential device applications in optics, 
microscopy, nanolith^raphjBSHSSl These studi es in the area of plasmonics also extend to carbon-based structures 
such as fullereneJi^^^, carbon nanotubeJ^^^^, as well as the recently discovered silicon-based silicene and other 
buckled honeycomb lattice structure d^^ l ^^ l with interesting roles played by the on-site potential differences between 
the sublattices. The main feature of the silicene energy dispersions is the energy gap, determined by the spin-orbit 
interaction. We believe that our theoretical formalism, here applied to gapped graphene, could also applicable in a 
straightforward way to silicene and germanene. 

One of our goals is to investigate the way in which the energy band gap modifies the temperature-induced plasmon 
excitations and their damping rates in intrinsic (or undoped) graphene. Specifically, we distinguish the cases when 
the graphene layer is free-standing and when it is in close proximity with a thick conducting substrate. Clearly, 
such plasmon modes cannot be excited at T = 0 because of the absence of free carrierJiSl The gaplessness of 
graphene is an important factor in the thermal population of the valence and conduction bands with electrons and 
hole correspondingly. Since the finite temperature also results in the decay into particle-hole pairs, we are also 
interested to study how this decay (damping rate) depends on the temperature. 

Closed-form analytic expressions for the long wavelength plasmon dispersion relation in gapped graphene as well as 
the dynamical polarization function were at T = 0 K initially reported by Pyatkovskiy in Ref. [S] . This important study 
demonstrated the existence of an extended range of wave vectors for the undapmed plasmons in the presence of a finite 
ener gy g ap. This bandgap could be opened using a substrate or by exposing the graphene layer to circularly-polarized 
light. ^ The polarization function in this case was obtained analytically in Ref. |22) . 

We also consider finite-temperature plasmons in graphene which is Coulomb-coupled to a semi-infinite conducting 
substrate. For possible tailoring of the plasmon frequencies, graphene has been combined with prefabricated plasmonic 
nanoarrays and metamaterials in order to obtain hybrid plasmon device^^lHlIl Therefore, a thorough understanding 
of the dispersion and damping of plasmons in graphene interfacing with different kinds of substrates is necessary 
for producing innovative practical applications. We have investigated zero-temperature non-local plasmons in one of 
our previous studie^^, applying the theore tical formalism for a graphene layer interacting with a surface plasmon in 
a semi-infinite conductor was developed irP^ESl. The plasma instability in such systems was addressed in Ref. [29]. 
These results could be considered as a non-trivial extension of Ref. [31] in which a linear acoustic plasmon mode was 
obtained for two interacting graphene layers. 
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Some of our work is devoted to a careful calculation of the analytical results for the real and imaginary parts of the 
polarization function in the long wavelength limit and for the dispersion equation which yields the plasmon modes. 
The rest of the paper is organized as follows: in Sec|^ we show that the finite-temperature polarization function 
used in Ref.[Tn] may be extended to the case of gapped graphene. After that, in Secjl^ we employ these results 
to gapped graphene to obtain the high-temperature plasmon dispersion relation in the long-wave limit, as well as 
the corresponding damping rate. This is an important contribution of our study. In Sec |III| we derive the finite- 
temperature plasma dispersion relation for a graphene layer which is Coulomb-coupled to a semi-infinite substrate. 
These analytical results show novel non-trivial behavior, connected with the temperature dependence of each plasmon 
dispersion for both real and imaginary parts of the plasmon frequency. Finally, numerical results for the plasmon 
dispersion using the full polarization function for gapped graphene are presented in Sec|IV[ 


II. HIGH-TEMPERATURE PLASMON DISPERSION RELATION FOR GAPPED GRAPHENE 

We now consider gapped graphene with energy dispersion: e(fc) = -|- (hvpk)'^ where vp the fermi velocity 

and A is the energy gap between the valence (-) and conduction (-I-) bands and the polarization function in the long 
wavelength limit is given by 
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Equation 0 is only valid for A < fi. If A > /r, the valence band is completely empty and only inter-band transitions 
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which is similar to Eq. (B5|. 


We now turn our attention to evaluate the finite-temperature polarization 
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which are still to be evaluated. 


A. Plasmon Dispersion relation (Re ojp) 
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which consists of two terms, namely, 
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Therefore, for T = 0 K, the first correction which is associated with the energy gap is linear, real and negative. We 
now collect all the approximations from each term and write 
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We note that the zero-order term is purely imaginary. We present the next part as 
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The remaining term is defined by 
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where G 2 (T, A) is a dimensionless integral defined as 
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The finite-temperature polarization function may be expressed as 
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From this, we deduce the temperature-induced plasma frequency of gapped graphene given by 
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B. Evaluation of Integral G2 

Here, we are looking for an approximated result in analytic form for the following integral: 
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where S = A/ksT. In order to avoid the singularity appearing when (5 —>■ 0, we perform the integration by parts. 
Consequently, we obtain 
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We estimate each term when d —> 0 or A <C ksT. We have 
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where Agk — 1-2824 is the Glaisher-Kinkelin constant, defined as 
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and 'jEu = 0.5772 is the Euler-Mascheroni constant, given by the following expression: 
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Numerical integration gives 
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and the damping rate ( 6 ) in units of 70 = rgV%k% log 2 ') for high-temperature plasmons in gapped graphene. 


Plasmon dispersions are shown for fcsT = l.OE^ andl.SE^'^ and energy gap values A = O.OandO.SE^ , and damping rates 
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and the plasmon frequency in Eq. (16) is now given by the equation 
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C. Plasmon Doping at a low temperature for extrinsic graphene 

Let us consider the case of a finite doping /i > 0 and low temperature T —> 0. We must calculate: the polarization 
function using 
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FIG. 2: (Color online) Finite-temperatnre plasmon excitations in graphene for various chosen values of doping and temperatnre 
(density plot of RPA polarization fnnction whose peaks correspond to the plasmon mode freqnencies). Panel (a) shows the 
zero-temperatnre limit of a graphene plasmon with finite doping /j, = 1.0E^\ Plot (6) corresponds to = l.OE^^ and 
H = l.OE^^ Panels (c) and (d) show plots of the plasmon dispersion relation for either high temperature or doping valne 
(fci,T = l.OE^^ and /r = 5.0E^^ for (c) and k^T = 3.0E^^ and fi = l.OE^' for (d). 


In the limit T —0, it follows that tanh(a/T) —sign(a), so when Eq. (27) gives 2 when A > /i and 0 when A < /r. 
So that the integral is 6{A — /i^. The second term is an integral representation of a Dirac delta-function for T = 0 K 
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D. Damping rates for gapped graphene 


The imaginary part of the zero-temperature polarization function is from Eq. 
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FIG. 3: (Color online) Overdamped (suppressed) plasmon frequencies for graphene at low doping concentration and T —>■ 0. 
Panel (a) gives the dispersion relation at low temperature fci,T = 1.0E^\ while panel (6) corresponds to T = 0. Both panels 
are plotted for intrinsic monolayer graphene. 
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This shows that when T = 0 K, the first-order correction associated with the energy gap is linear, real and negative. 
At finite temperature, our calculation shows that 
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Additionally, making use of similar procedures when we dealt with zero gap, we obtain 
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from which we deduce the damping rate 
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FIG. 4: (Color online) Density plot of the imaginary part of the polarization function at various chosen temperatures and 
doping concentrations in graphene. Panel (a) shows a plot corresponding to low temperature kpT — O.OIE® and zero doping 
(in which plasmons do not exist due to the absence of free carriers). Plot (b) demonstrates the low-temperature limit for finite 
fM = 1.0E^\ Panel (c) corresponds to finite temperature and zero doping, whereas panel (d) shows the finite-temperature 
particle-hole modes for highly doped graphene n = 5.0E^\ 


III. GAPLESS MONOLAYER GRAPHENE INTERACTING WITH A THICK CONDUCTING 

SUBSTRATE 

According to recent worlspfll23j the plasma dispersion relations when a 2D layer is Coulomb-coupled to a thick 
conducting substrate are determined by solving the following equation: 
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where a is the layer-surface separation, ujp is the bulk plasmon frequency and the Fourier transform of the Coulomb 
potential energy 
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is expressed in terms of Vg is a dimensionless parameter. 
The finite-temperature polarization function is 


27r 


hrgVp 


q 


(40) 


4°^(g,u;) 


21 n2 

TT 


kpT + 


16 keT 


(41) 


Here, we consider a well-defined plasmon with j uj, which is excited with a certain wave vector at a chosen 
temperature. Accordingly, we neglect the contributions to the real part of the polarization function arising from the 
imaginary part of the frequency. So, making use of 
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FIG. 5: (Color online) Plasmon dispersion relations for monolayer graphene, which is Coulomb-coupled with a semi-infinite 
conductor for different chosen values of graphene doping and temperature. Panel (a) shows the low-temperature limit when 
the doping is finite corresponding to /r = 1.0E^\ Plot (&) and shows the case for fci,T = O.IE^^ for a slightly doped graphene 
^ = 0.1E^\ Both panels (c) and (d) correspond to a moderately high temperature ksT = 3.0E^^ and different chemical 
potentials {fj. = l.OE^^ for (c) and /i = 5.0E^^ for (d)). 


in Eq. (39), we obtain 
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We now consider two different limiting cases. The first corresponds to when the surface-layer separation is small so 
that aq <?; 1. For this, we have the results in the linear approximation 
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Both solutions are linear in q, we also note that the surface plasmon frequency does not depend on a until the second- 
order terms for the 17 1 . branch. Second-order terms in q are neglected since we are interested in the temperature 
dependence in the lowest order. Turning now to the set-up when aq^ 1, a straightforward calculation yields 


In the local limit Ag ^ 1, we obtain 


111 

UJp 

Q 2 

UJ. 




P 


21-Xq 
1 Xq 


■\/2 \/2 1 — Xq 


(46) 






172 = -^+2^/2\a2^^kBT q& 

V2 Aicjp 


-2aq 


(47) 


We need to take into account the imaginary part of the — ojp) term. We make the replacement uj ^ oj + ij 

which leads to the following plasmon dispersion equation 
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When the surface-layer separation is large, the imaginary parts may be further approximated as 


(50) 


hi 

r 2 








l-81n2^^fcBrqe-2“9 


huj 


t -2aq 

- huJr> - QQ ^ 

16 


(51) 


where A has been dehned in Eq. (44). We now investigate how some of these results are affected when a finite energy 
gap exists between the valence and conduction bands. 


IV. PLASMONS IN GAPPED GRAPHENE, COULOMB-COUPLED TO A SEMI-INFINITE 

CONDUCTOR 

The real part of the polarization function may be expressed as 
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on the other hand, when aq 3> 1, the solutions are given simpler are given by 
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We now turn to presenting and discussing our numerical results in the next section. 
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V. NUMERICAL RESULTS AND DISCUSSION 

The plasmon dispersion relations and damping rates for gapped graphene in the long-wavelength limit are presented 
in Fig. The plasmon frequency is increased at higher temperature, approximately having the behavior of 
77 dependence as for gapless graphene. Clearly, the damping rates decrease with increasing energy gap for all 
temperatures. This means that the plasmon modes become less damped in the presence of a finite energy gap. 

We are definitely interested in calculating the polarization function and plasmon energies numerically beyond the 
long wavelength limit. First, we obtain the imaginary part of the polarization function in Fig. looking for the 
regions within {q — wj-plane where the relatively undamped plasmon modes exist (im n*^°)(q,a;) —>-0). As we know 
for the long wavelength limit, the imaginary part of the polarization function within the region of graphene plasmon 
excitations decreases at high temperature as 1 so it is relatively small at high temperatures although it not be 
negligible. Panel (&) of Fig. ^demonstrates the behavior when/r = at very low temperature, which approximately 
corresponds to the standarcTparticle-hole modes of zero-temperature graphene.!^ We see that at finite doping and high 
temperature, the region of particle-hole modes (finite Im (q, w)) above the diagonal w = vpq is suppressed, as 
shown in panel (d). 

Figure [^presents the RPA polarizability for various doping concentrations and temperatures. At zero temperature 
and when the doping is finite, we reproduce the well-known plasmon dispersion relation in graphene, as it was reported 
in Ref. [ 3 . As the temperature is increased, the plasmon frequency is also increased for all values of the wave vector, 
growing with temperature like 77 in the long wavelength limit as it was showed analytically in Ref. m)- Various 
examples of damped plasmon excitations are shown in Fig. As it was mentioned above as well as in Ref. cni, the 
plasmon mode has low intensity in the limit of vanishing temperature and zero doping. 

The temperature-dependent plasmon dispersion relation for a graphene monolayer interacting with a semi-infinite 
conductor is shown in Fig. There are clearly two plasmon branches, originating from the graphene layer (acoustic 
branch, starting at the origin), and the “sruface plasmon”, which is depolarization shifted from its long wavelength 
value of at l/TSwp. Both branches have positive group velocity for {q —> 0), i.e. are linear in the long-wave limit, 
similarly to its variation at zero temperature. Another interesting detail, which deserves mentioning is that the higher 
plasmon branch, attributed to the surface behaves similarly to the plasmons in graphene (see Fig.[^, while the lower 
acoustic branch remains nearly unchanged for all the chosen temperature and doping values. 

The plasmon frequencies and the damping rates for grpahene. Coulomb-coupled with a thick conductor with free 
carriers in the long wavelength limit are presented in Fig. Both plasmon dispersions are linear and increase with 
temperature as ire 77 and ire T, respectively. We also note that the damping rate F 2 , which corresponds to the upper 
plasmon, depends linearly on the wave vector q. This is a new result, which was not encountered in graphene, either 
with a finite or zero energy gap. 
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FIG. 6: (Color online) (a) Plasmon frequencies and (6) the damping rates in the high-temperature limit for monolayer graphene 
interacting with a semi-infinite conductor in the long wavelength limit. Each plasmon branch (both real and imaginary parts) 
are presented at two chosen temperatures - ksT = l.OE^^ and fesT = 2.0E^\ 


VI. CONCLUDING REMARKS 

In summary, we have obtained analytic expressions for the high-temperature plasmon dispersion relations for gapped 
graphene, as well as for a graphene layer interacting with a semi-infinite conductor. We have found that the plasmon 
frequency is modified according to /{AksT) ln(A/ {2kBT)) in the presence of an energy gap, which is different from 
the case when T = 0. The corresponding damping rate is decreased linearly in the presence of a gap, making the 
plasmon more stable compared to its counterpart for gapless graphene. Our investigation of a graphene monolayer, 
which is Coulomb-coupled to a semi-infinite metal, revealed a number of novel features for both plasmon frequencies 
and the damping rates. Specifically, we emphasize the linear (/—dependence of some of the damping rates, which has 
not been encountered for graphene (either gapped or ungapped). 
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Appendix A: Finite temperature polarization - integral transformation derivation 


In this Appendix, we pr ove that the finite-temperature polarizability could be expressed as an integral transfor¬ 


mation given in Eq. 

Ref. [To]. 


(All relating the polarization function n^^((/,w) at zero temperature, as it was reported in 
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We need to prove the general validity of the following expression: 







AkeT 


^cosh 


2kBT 


2 


with the zero-temperature polarization function is as follows: 


(Al) 




np[e'^{k)] - np e'’'(|k + q|) 


e®(fc) - e®'(|k + q|) -|- h{uj + * 7 ) 


(A2) 


in terms of th e sp in and valley degeneracy factor g = A and the form factor (q, k) . The only temperature-dependent 
terms in Eq. (A2) are the Fermi-Dirac distribution functions ^^(e) = (l-|-exp[(£— p,)/(A:bT)])“^ = (l-t-exp[/?(e—^)])“^ 
with /3 = l/^ksT). This means that we must prove a similar integral transformation for the distribution functions. 
The Fermi-Dirac distribution functions could be presented as 


np{0 


1 

2 


— tanh 



(A3) 


where ^ = e — g. Indeed, 

2 - 1 _ g/3«/2 

^ ~ 1 -Pe« ^ e/3? -p 1 ^ e/3?/2 e/3?/2 ^ T 


(A4) 


While at zero temperature, the distribution function could be presented as the Heaviside unit step with np{e — g,,T —>■ 
0) = 9{p. — e). Accordingly, 


00 


0 


9{fi — e)dg' 

AkpT cosh^lifi - g')/{2kBT)] 


00 



€ 


d-d' 

2kBT 


(A5) 


We substitute rj = 2 k bT dg! = —2kBTd'q. The limits of integration now become ry< = 2 kB^T 9> “ 

Obviously, ??< > rjy. The integral now takes the form 


Also we can write 


— 00 


/ 


—2kBTdri 
AkBT cosh^ [ij] 


v< 



— 00 


1 

2 


(tanh77< -|- 1) . 


Finally, we obtain 


tanh r]^ 


tanh 


/i — e 
2kBT 


— tanh 


2 


{e- g) 


np{0 


1 

2 


1 - tanh ^ 
2 


Alternatively, we have 



_ did - g) _ 

AkBTcosh^lifi - ti')/{2kBT)] 


00 

(p ■ T\= f np{£, g'\T = Q) dg' 

7 4A:BTcosh2[(/r-^0/(2fcBT)] ■ 


(A6) 


(A7) 


(AS) 


(A9) 






























Appendix B: Derivation of the long wavelength limit of the Polarization FunctiorP® and the 

Finite-temperature Plasmon Moded^ 


First we need to derive the long-wave limit of the polarization with g —>■ 0 with hxed frequency uj vpq. We 
consider two separate cases with different imaginary part of the polarization: 1) w > 2/r and 2) a; < 2/r. 


n°(g,cc) = 


^ , 1 In 

SttHu! \ hu! 2 


2/r — Huj 


2fi — huj 

If we consider the region with huj <C ft, we can simplify as follows: 


- y - 2ft) 


In 


2fi — fio; 


Huj 

2ft 


2fi — Huj 

Finally, for the region lA (where the zero-temperature plasmons exist), we have 


gg^ 2ft ^ ft 
SttHlo Huj ttH'^ uj'^ 

Im n°(g,a;) = 0, {Hw < 2ft) 


Re n°(( 7 ,w) = 


(Bl) 


(B2) 


(B3) 


1. Derivation of Eq.( Bl) 


The general exression for the non-interacting polarization functions is: 

n°(g,a;) = Po{q,u}) +SP{q,Uj) 


(B4) 


with the term Po(q,w), corresponding to the inter-band transitions, which is urely imaginary for oj > vpq and the 
fi-dependent 5P{q,uj)^ which appears due to the inrta-band transitions inside the conduction band (for ft > 0). The 
terms could be introduced in the following form 


Po(g,w) = 

6P{q,uj) = - 


(B5) 


Jp 

gM 


2'Kh‘^v'jp h'^vp 


1 [ fHuj + 2fi\ /2ft-/tw 

Fi{q,uj) { Fa ( ) - -^2 I -A - I - 


Hvpq 


Hvpq 


with 


hvW 


^i(g,w) = — -^ 

lOTT 

F 2 {X) = X^X‘^ - 1 - ln(A’ -b - 1) (T” > 1) 


(B6) 


and Vpq <C w. We start with the following approximation 


Eq.(B5) is satisfied in the Region lA, given as w < 2ft — vpq and uj > vpq. Now let us consider each term for q ^ 0 

(B7) 


1 




- 1/2 


1 (vpq)^ Svj.q'^ 


yjuj’^ — v^q^ w \ \ UJ J J UJ 

Consequently, the inter-band polarization function has the form 

mq^ ip v\q'^ 


2a;3 


8 uj^ 


(B8) 


i.e. is purely imaginary. 

Now we consider the intra-band part 6P{q,uj) . First, we analyze how each term of it behaves for q —>■ 0. Thus 
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and 


F2{X) = X\JX’^ - 1 - ln[A' + x/A’^ - 1], 

A’a/A’2 - 1 = A’2(1 - IjX'^Yl^ - A’^ - i 

A- - VA’2-1 = A-^ - a- ("l - ^X^-X 


Let us summarize the results: 

i' \ / \ / \2 / \2 

hu! + 2^\ j huj — 2fi\ jHu! + 2^\ 1 j hu! — 2iJ, 


' hvpq 
2 X Ahojn 
{hvpqY 

and also 


As a result 


-Fn 


2 In 


hvpq 

2fi — hoj 
2 /i + hui 


\ hvpq 


2 \ hvpq 


+ --ln 


^ T? i ^ 9 q 

Fi{q,uj) x: 




2 4 

_ 9 vpq 

16ttH oj 327rS w 2 




Fi{q,uj) 


Hij A 2^ 

hvpq 


-Fo 


/ So; — 2/i 
\ hvpq 


99 


(A> 1) 


/ fiw + 2/i 
\ hvpq 


+ ln 


_ 99 9 

2TrlFvjp 4fi,2 7ra;2 


(BIO) 


/ huj — 2/i 
\ hvpq 


(Bll) 


(B12) 


(B13) 


Finally, we can write 



1 In 

2fj, — huj 

\huj) 

Idirhu; 

2ijl + huj 


xuf \ 99 , 99 , 9 

As far as the imaginary part is concerned, it is equal to zero everywhere in Region lA: 

Im dF[q,uj) = nr — 


According to Eq. (B5) 


so that 


ITT 

Po{q,u}) = -^^Fi{q,uj) 


n^{q,u;)=P°{q,Uj) + 6P{q,Uj)=0 


(B14) 


(B15) 


(B16) 


(B17) 


So Eq.(Bl) is confirmed. 


2. Finite-temperature plasmons 

The finite temperature non-interaction polarization function is expressed as follows: 

OO 

n°(g,a;) 


n®(q,a;) = J ■ 


2kBT 


iksT ^cosh 

In the long-wave limit approximation, the zero-temperature polarizability is as follows: 


rrO^ ^ 99^ 9 ■ 99^ n/t o ^ 

n (q,^) = - ^9) 


47r (Scli)2 Idhuj 


(B18) 


(B19) 
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Let us first consider the real part: 


Re = f _^_ ^ _ 

^ 7 47r (fia;)2 AkeT cosh^ [{fi - ^i')/{2kBT)] 


This integral could be easily evaluated: 


gq^ 1 

Airh^u)'^ AksT 


oo 


djj! 


cosh^ 


2kBT 


(B20) 


dfi' 


cosh^ 


2kBT 


= Ak'^T'^ In [1 + 


We consider intrinsic graphene with /i = 0 so that In [l + = In 2. 

Thus 


T5 tt(o)/ ^ 5ln2 

Re = 


The imaginary part is as following: 


(B21) 


(B22) 


Im n<?>(,.,.) = I 


, gq^ f dg' 6{huj — 2g') 


eosh^ 


M-M 


2kBT 


(B23) 


Finally 




Whoj 2 AkuT 128 kuT 


(B24) 


3. plasmons 


Now we rewrite the finite-temperature polarization function: 


( 0 )^ 21n2 ^ ^ ^ i 

(,, 0 .) = —^t«r+- — 


The dielectric function e{q, lo) is 


(B25) 


e(g,w) = 1 -V{q)U^°\q,uj) 


(B26) 


where V{q) = 27re^/(esg). 

We introduce real and imaginary parts of the frequency w —>■ w -I- iy and /{eshvp), we find the real part: 


2716^2 In 2 q'^ , rj. a 
1 - ^^^bT = 0 

Csq IT 


(B27) 


we obtain: 


2 41n2 

w = , vprsq ksT 


(B28) 


Appendix C: Imaginary part of the frequency 

For A = 0 the polarization function is 




21 n2 q^ i q^ 

TT s + 16 keT 


(Cl) 



The plasmons are defined by the following equation 
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e{q,ijj) = 1-n^^(g,CLi) = 1- rshvFll^\q,u;) = 0 

Q 


Let us introduced the complex frequency: 


w uj + i"f 


up {uj + i'^Y 


1 - 2 ^ 7 / 0 ; ^ _ 2 ^ 7 


This results in the following equation 


27r 


According to Eq.(B28) 


which leads to 


and 



2^7^ 


_ 

a;3 J 

' 16 keT j 


2 In 2 

Tr/i^ 


27 r . f 2 In 2 9, 1 . 

1 - TshvF { )■ = 0 

q I irn^ ui^ 


= 0 


2 In 2 2 , rrPh , * 9^ 




q ksT —F + 


7 = 


a;3 lefcsT 


ttIP a;3 


= 0 


In 2 64yfc2r2 

^3 = 


TT ! ^3/2 

7 = ^ (VFrsq) 
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(C2) 

(C3) 

(C4) 

(C5) 

(C6) 

(C7) 

(C8) 

(C9) 
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